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Developments in the evaluation of rric dependent matrix elements in 
B -^ X,7 at NNLO 
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Am Hubland, D-97074 Wiirzburg, Germany 

The full prediction of the inclusive radiative B — > Xs-f decay rate at NNLO requires a complete evaluation of 
missing charm quark mass dependent matrix elements. Recent developments in the on-going computations are 
reported and the current status is briefly overviewed. 



1. Introduction 

One interesting candidate process in the in- 
direct search for non-standard physics is the 
rare inclusive B -^ Xgj decay [1]. Due to 
its low sensitivity to non-perturbative effects 
and, being a flavour-changing neutral current, 
its loop-suppression in the Standard Model (SM) 
stringent constraints on the parameter space of 
physics beyond the SM can be derived from both 
accurate measurements and precise theory pre- 
dictions. 

The latest measurements by BaBar, Belle and 
CLEO [2] have been combined by the Heavy 
Flavour Averaging Group (HFAG) [3] into the 
current world average (WA) for the branching ra- 
tio. 



„exp 



(3.15 ±0.23) X 10" 



(1) 



where a cut E^^q > 1.6 GeV has been imposed 
on the photon energy in the i?-meson rest frame. 
The first uncertainty corresponds to a combined 
statistical and systematical error, the second one 
is due to the theory input in the extrapolation 
of the measured branching ratio to the reference 
value Ej^o, whereas the third one is connected to 
the subtraction of 6 — > ^7 contamination. On the 
theory side, the expected size of next-to-next-to 
leading order (NNLO) QCD effects to the par- 
tonic decay b -^ ^partomc^ jg comparable with 
the overall experimental error of about 7% and 
therefore a complete SM calculation at this level 



of accuracy is clearly needed. 

The recent theoretical estimate of the branch- 
ing ratio at the NNLO level 



^f;>i.9GeV = (3.15 ± 0.23) X 10-4 



(2) 



was derived in [4J after a large part of the NNLO 
program has been finished and is in good agree- 
ment with the WA. The uncertainty here consists 
of four types of error added in quadrature: non- 
perturbative (5%), parametric (3%), higher-order 
(3%) and TOc-interpolation ambiguity (3%). 

Large logarithms of the form 

as{mb)'^\o^{mb/mw) appear in QCD correc- 
tions to the partonic decay width r(6 — > 57) and 
have to be resummed with renormalization-group 
techniques to get a reasonable prediction. Most 
suitably, this is done in the framework of an ef- 
fective low-energy theory with five active quarks 
by integrating out the heavy electroweak and the 
top fields in the SM. As a consequence, local 
flavour-changing operators Qi{fJ-) up to dimen- 
sion six and Wilson coefhcients Ci{ii) appear in 
the resulting effective Lagrangian £eff ■ 

A few years ago, the next-to-leading order QCD 
corrections to the 6 — > 57 decay have been com- 
pleted (see e.g. |5I6) and references therein). The 
next-to-next-to leading order evaluation, which is 
a very complicated task, is currently under way 
and large parts are already finished. In general, 
three steps are required for a consistent calcu- 
lation in the low-energy effective theory and in 
particular also at the NNLO level: 



1. Determination of Wilson coefficients Ci{fio) 
at the electroweak scale fiQ — Mw by re- 
quiring equality of Green's functions in the 
effective and full theory at leading order 
in (external momenta) /Mvi^. To this pre- 
cision, the matching of the four quark op- 
erators Q17 ■ ■ ■ 7 Qe and the dipole operators 
Qy, Qs at the two- and three- loop level, re- 
spectively, has been computed in |7l8j . 

2. Derivation of the effective theory Renor- 
malization Group Equations (RGE) and 
computation of the operator mixing under 
renormalization by evolving the Wilson co- 
efficients Ci(/Lt) from Ho down to the low 
scale fXb ~ mb using the anomalous dimen- 
sion matrix up to 0{al). In the sectors 
{Qi, . . . , Qe} and {Qy, Qs}, the three-loop 
renormalization was found in |9I10| . Results 
for the four-loop mixing of {Qi, . . . ,Qe} 
into {Q7, Qs} were recently found by com- 
pleting the anomalous dimension matrix 

3. Calculation of on-shell matrix elements 
with single insertions of effective operators 
at /if, ~ rub to 0{al). This task is not 
completed yet, although a number of con- 
tributions is known. In [12|13j . the two- 
loop matrix element of the dipole oper- 
ator Qr together with the corresponding 
bremsstrahlung was determined, confirmed 
in [14] and subsequently extended to in- 
clude the full charm quark mass depen- 
dence in [15j . Dominant contributions to 
the photon energy spectrum in the so-called 
large-/3o approximation 0(1300^) have been 
obtained in [16j . Using an expansion in 
ml/ml the 0{Poa1) contributions to the 
two- loop matrix elements of Qj and Qg) as 
well as to the three-loop matrix elements 
Qi and (52, were found in [17]. These re- 
sults have been confirmed in [18j and, more- 
over, the full fermionic corrections beyond 
the large-/3o approximation have been pro- 
vided there. 

The full matrix elements of Qi and Q2 at 0{a1) 
constitute an important piece that is still missing. 



At NLO the choice of scale and scheme of TOc 
constitutes the main source of uncertainty stem- 
ming from the fact, that these operators contain 
the charm quark and contribute for the first time 
at 0(0;^). Removing this ambiguity is therefore 
a NNLO effect in the branching ratio. So far, 
the full matrix elements of the Qi and Q2 oper- 
ators have been evaluated in the large rric limit, 
TOc 3> mb, and subsequently used for an inter- 
polation to the physical range of TOc [TSj assum- 
ing some ad-hoc value at rric = 0. This is the 
source of the aforementioned interpolation ambi- 
guity in the current NNLO branching ratio esti- 
mate. Removing this uncertainty requires the cal- 
culation oi {s^\Qi^2\b) at physical TOc and involves 
the evaluation of hundreds of on-shell three-loop 
vertex diagrams with two scales rrib and TOc- Re- 
ducing this interpolation uncertainty is possible 
by computing the matrix elements at rric = 
and thus fixing the endpoint of the interpola- 
tion. Both evaluations are currently being pur- 
sued. This paper is meant to describe the current 
status. 



2. The matrix elements {s^\Qi,2\b) at 
NNLO 

The effective low-energy theory, in which the 
calculation of {s^\Qi^2\b) at 0{a1) is performed, 
is given by the effective Lagrangian 

CcS — ^QCDxQEDiu, d,S,C,b) 

+^F,:y,,^a(A*)g.(A.). (3) 

^ i—l 

The first term corresponds to the usual QED- 
QCD Lagrangian for the light SM fields, the sec- 
ond term gives the local operator product expan- 
sion, Vij are the matrix elements of the Cabibbo- 
Kobayashi-Maskawa matrix and Gp is the Fermi 
coupling constant. The relevant physical opera- 
tors are chosen as [2S1 



Qi,2 = (sr,c)(cr^6), 

(53,4,5,6 = (sri6)Y] (gE-q), 



(4) 



V7 = 

Qs = 



■rnbi^i)isL(J^''bR)F^, 



p.i'j 



167r2 



where F and V stand for various products of 
Dirac and colour matrices. 

As already mentioned, removing the interpola- 
tion uncertainty requires the calculation of Qi^2 
matrix elements at the NNLO level. All ap- 
pearing vertex diagrams have been generated, ex- 
pressed through scalar diagrams that depend on 
the two mass scales rric and rrib and finally re- 
duced with Laporta's algorithm to 476 master in- 
tegrals. The latter can be evaluated using two dif- 
ferent approaches: numerical solutions of differ- 
ential equations or the Mellin-Barnes technique. 
The first method consists of writing down differ- 
ential equations for the master integrals in their 
kinematical invariants and solving them numer- 
ically after an expansion about some limit has 
been taken to provide the startpoint for inte- 
gration with high precision. This was already 
used in the calculation of double fermionic QCD 
corrections to the photon polarization function 
|20j . Moreover, in |21j this idea was also ap- 
plied to determine the full mass dependence of 
the it-production cross-section from light quarks 
at NNLO. The second approach is based on the 
Mellin-Barnes (MB) technique, where MB repre- 
sentations can be derived in an automatized way, 
analytically continued in e and numerically inte- 
grated utilizing the package MB [52]. Both meth- 
ods have been used to obtain the full fermionic 
corrections, extending the so far known massless 
case to include heavy b and c quark loops [TB] . As 
a result, it turned out that while the charm quark 
contribution is well reproduced by the massless 
approximation, bottom quark loop insertions are 
overestimated by a large extent. For the remain- 
ing bosonic parts the IBP reduction is still under 
way. In the limit rric — valuable information 
can be obtained from {sj\Qi_2\b) to fix the end- 
point of the interpolation performed in [19|. A 
possible way of getting {sj\Qi^2\b) is by interfer- 
ing the matrix elements of the operators (5i,2 with 
those of the dipole operator Qy and cutting the 
resulting four-loop propagators selectively. The 
final state is required to contain at least one s 




Figure 1. Example graph 

quark and one photon. Altogether, several hun- 
dreds of four-loop propagator-like diagrams are 
generated with up to five-particle cuts and sub- 
sequently reduced to about 200 master integrals 
that have to be evaluated. The massless cases 
among them have already been computed up to 
four-particle cuts using a Mellin-Barnes based 
method. For this task, Mellin-Barnes represen- 
tations have been derived and numerically evalu- 
ated after performing the necessary phase space 
integrations. However, integrals containing mas- 
sive lines from virtual 6-quarks require a different 
approach. The method we are currently using 
in this respect is again based on differential equa- 
tions (DEQ), but it turns out to be more involved 
as compared to the virtual corrections at rric ^ 0. 
As starting point for the expansion and integra- 
tion of the system of DEQ we choose the large 
mass limit. Thus, boundary conditions are given 
by cut integrals in the large mass limit and can be 
evaluated with automatized diagrammatic expan- 
sions and phase space integration. Unfortunately, 
the usual approach of numerical integration up 
to the interesting kinematical point, to the on- 
shell limit in our case, suffers from divergences 
in the differential equations at both endpoints of 
the integration contour. This can be overcome 
by solving the system of DEQ with an expansion 
ansatz around the on-shell point. Comparison of 
the numerical integration at some point close to 
the on-shell kinematics with the expansion fixes 
the boundary conditions of the latter. In the case 
of logarithmic on-shell divergences, the integrals 
are regularized with a sufficient number of irre- 
ducible numerators and a change of the basis of 
master integrals. 

The achievable accuracy of this method is, on 



the one hand, given by the depth of the expan- 
sions at the starting and ending points on the in- 
tegration contour and, on the other, by the error 
pile-up during numerical integration. Performing 





: p^ /m? 



Figure 2. Example integral together with its plot 
as function of the ratio z = p^ /m?. Dots de- 
note the exact numerical results and the lines at 
z < 0.5 correspond to the series expansion for 
z ^ 0. Curves at z > 0.5 show different depths 
in the expansion around the on-shell point, the 
matching is performed in z = 0.9 

the numerical part with quadruple precision, a 
relative error of the order of 10^^^ at the on-shell 
point seems to be obtainable. 

As a first attempt we study diagrams involv- 
ing two- and three-particle cuts, which constitute 
the major part of the evaluation. Fig. [2] shows 
a plot of the ©(e") term for one example inte- 
gral. The expansion in the large mass limit has 
been performed up to 20 terms in z = p^/rn^ 
whereas in the on-shell limit 12 terms have been 
taken into account. The numerical integration is 
performed starting at z = 0.05 and the match- 
ing of the numerics and the on-shell expansion is 
done dX z — 0.9. It is apparent, that both series 



converge nicely against the numerical points. 

The first three terms of the expansion in the 
limit oi z — p^ /m? — > of the integral in Fig. [2] 
are given by 



-fz^O 



10 

T 



2\og{z) 



3 V 4 2 2 ' 



h -C2 log(2;) 

54 9 3 



98 



+ 4C2 



28 

3 ■ '" "g" 

31og(z)3 + z 



C3-341og(z)-151og(2)2 



51 21^ 5^ 
y - ^C2 + -C3+ 



(-f -^C2)log(.)-^Og{.)2J + 

2047 13 ^ 10 ^ 

-^ + |C2Jlog(^)-ilog(2)2 



The corresponding first few terms of the expan- 
sion around the on-shell kinematical point with 
y = z^^ — 1 read 

h^i = 



-2.0944- 



■ (-12.703-1-10.335?;- 



4.7124 y'^) + - (-52.607 + 81.505 y - 67.338 y'^) 
+ (-175.32 + 454.591/- 472.67 y^) . 



3. Conclusions 

An evaluation of the missing matrix elements 
{sj\Qi,2\b) at 0{af) is essential for the reduction 
of the current uncertainty in the estimate of the 
branching ratio B{B — > X^j). Two different ap- 
proaches are used and will eventually reduce or 
even remove the remaining ?7ic-interpolation am- 
biguity that amounts to about 3%. 
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